The mean-field y> 4 -model: entropy, analyticity, and configuration space topology 



Ingo Hahn* and Michael Kastncr^ 
Physikalisches Institut, Lehrstuhl fur Theoretische Physik I, 
Universitat Bayreuth, 95440 Bayreuth, Germany 

A large deviation technique is applied to the mean-field (/? 4 -model, providing an exact expression 
for the configurational entropy s(v, m) as a function of the potential energy v and the magnetization 
m. Although a continuous phase transition occurs at some critical energy v c , the entropy is found 
to be a real analytic function in both arguments, and it is only the maximization over m which gives 
rise to a nonanalyticity in s(v) = sup m s(v,m). This mechanism of nonanalyticity-generation by 
maximization over one variable of a real analytic function is restricted to systems with long-range 
interactions and has — for continuous phase transitions — the generic occurrence of classical critical 
exponents as an immediate consequence. Furthermore, this mechanism can provide an explanation 
why, contradictory to the so-called topological hypothesis, the phase transition in the mean-field (p 4 - 
model need not be accompanied by a topology change in the family of constant-energy submanifolds. 
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I. INTRODUCTION 

The theory of phase transitions is an active field of 
research already since many decades. Despite several re- 
markable achievements on both the conceptual and the 
applied side, there still remain open questions galore. 
As a contribution to the conceptual side, a recent and 
promising proposal is the topological approach to phase 
transitions, connecting the occurrence of a phase transi- 
tion to certain properties of the potential energy func- 
tion vn, resorting to topological concepts. Among the 
remarkable features of this approach we mention: 

1. Topology is in some sense a very reductional level 
of description. Compared to a geometric descrip- 
tion of the same object, lots of information is disre- 
garded in topology. So the topological approach is 
an attempt to trace back the occurrence of a phase 
transition to more fundamental quantities than the 
measures in state space commonly considered. 

2. The microscopic Hamiltonian dynamics of a system 
can be linked via the maximum Lyapunov expo- 
nent to topological quantities [1]. These topologi- 
cal quantities, in turn, arc linked to the occurrence 
of a phase transition by the topological approach. 
In this way, a connection is established between 
phase transitions on the one hand and the under- 
lying microscopic dynamics of the system on the 
other hand. 

The topological approach is based on the hypothesis [2] 
that phase transitions are related to topology changes of 
a family {M v } of certain submanifolds M v of the configu- 
ration space of the system. For a system with N degrees 
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of freedom, the M v consist of all points ip of the configu- 
ration space r for which the potential energy V]y((p) per 
degree of freedom is smaller than or equal to a certain 
level i>, 

M v = {ip€T\v N ((p)^v}. (1) 

(Or, in a related version, the topology of submanifolds 
£„ = {ip £ T | vn({p) = u} is considered.) The hypoth- 
esis then conjectures that a topology change within the 
family {M v } at v = v c is a necessary condition for a ther- 
modynamic phase transition to take place at v c [or at the 
corresponding critical temperature T c = T(v c )]. This hy- 
pothesis has been corroborated by numerical and exact 
results for a model showing a discontinuous phase tran- 
sition [3, 4] as well as for systems with continuous phase 
transitions [1, 5-11]. A major achievement in the field is 
the recent proof of a theorem, stating, loosely speaking, 
that, for systems described by smooth, finite-range, and 
confining potentials, a topology change of the submani- 
folds M v is a necessary criterion for a phase transition to 
take place [12, 13]. 

Recent results on the configuration space topology of 
the mean-field ip 4 -modc\ have cast some doubt on the 
general relevance of topology changes for the occurrence 
of phase transitions. It was first observed in Ref. [14] 
that, for this model, the potential energy v c at which the 
thermodynamic phase transition takes place does not co- 
incide in general with the potential energy corresponding 
to any of the topology changes of {M v }. Independently, 
this result was confirmed in Ref. [15]. The discrepancy 
between the critical energies from thermodynamics and 
from topology provoked quite some speculations in the 
literature [14, 16]. 

In the present paper, an explanation for this discrep- 
ancy is given. The explanation is based on the observa- 
tion that a topology change in the family {M„} is one 
possible mechanism to entail a thermodynamic singu- 
larity. The theorem of Refs. [12, 13] then states that 
for systems with smooth, finite-range, and confining po- 
tentials, a topology change is the only mechanism to 
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generate such a singularity. We individuate a further 
singularity-generating mechanism, namely a maximiza- 
tion over one of the variables of a real analytic func- 
tion [31], which occurs in the mean-field (^ 4 -model. The 
availability of this mechanism accounts for the fact that 
the phase transition in this model need not be, and in 
fact is not, related to the topology changes observed. 
This second singularity-generating mechanism is by no 
means restricted to the mean-field < / 9 4 -model, but is gen- 
uine to systems with long-range interactions, impossible 
to take place in short-ranged ones. Not only can this 
maximization over a real analytic function account for 
the absence of a relation between thermodynamic and 
topological quantities, it also leads generically to the so- 
called mean- field (or classical) critical exponents a = 0, 
(3 = |, ... in case of a continuous phase transition. 

To the purpose of illustration of this singularity- 
generating mechanism, we present an analysis of the 
microcanonical entropy function of the mean-field <p - 
model, which could stand also as a result of interest on 
its own. The entropy s is derived from a large deviation 
principle, a powerful tool from probability theory. As a 
function of the potential energy v and the magnetization 
to, the entropy is found to be real analytic, and only the 
maximization over to generates a nonanalyticity in s(v). 

After fixing some notation in Sec. II, the microcanon- 
ical entropy s of the mean-field (/j 4 -model is computed in 
Sec. III. This result is confronted with the canonical one 
in Sec. IV, where in particular the (non-)equivalence of 
the statistical ensembles is investigated. In Sec. V, part 
of the results on the topology of the M v from Refs. [14-16] 
is reviewed, and the relation to thermodynamic quanti- 
ties is discussed. Section VI is devoted to two different 
ways in which a thermodynamic singularity can be gener- 
ated: by topology changes in the M v , or by maximization 
over one of the variables of a real analytic function. In 
the final Sec. VII, we comment on the general applicabil- 
ity of our results to long-range systems, as well as on the 
altered perspective from which the topological approach 
to phase transition has to viewed as a consequence. 



II. MEAN-FIELD V 4 -MODEL 

The Ginzburg-Landau or field-theoretic </? 4 -model 
plays a prominent role in the theory of critical phenom- 
ena, mainly due to the fact that it provides a unifying 
description of various, physically different, lattice and 
continuum systems near criticality (see [17] for an intro- 
duction) . Here we study a lattice version of this model, 
where the degree of freedom associated with the i-th. lat- 
tice site is characterized by a scalar variable tpi € 1R. 
As an analytically tractable caricature of physically rel- 
evant long-range forces, we consider mean field- like in- 
teractions, where all degrees of freedom interact with 
each other at equal strength. For such a system consist- 
ing of N degrees of freedom, the Hamiltonian function 



H N ■ H -> H is given by 
1 N 

H N (7r,<p) = -J2^ + N[v N (<p)}, (2) 

i=l 

where J is a coupling constant, and n = . . . , 7Tjv) is 
the vector of momenta conjugate to ip = (tpi, ■ ■ ■ , <Pn)- 
The potential energy is given by 

v N (<p) = z N (<p) - | [m N (<p)] 2 , (3) 

with on-site potential 

i=l v ' 

and magnetization 

1 N 

m N (ip) = —^cpi, (5) 
i=i 

where vn, zjv, and tun are functions from R w onto the 
reals. Typically, a coupling term —hm((p) to a symmetry 
breaking magnetic field h is added to the Hamiltonian 
(2). Since we arc interested in the zero- field situation, 
such a term will here be omitted. 

Canonical thermodynamic quantities of the mean-field 
<yj 4 -model have been computed in Refs. [18, 19]. We will 
make reference to these results when discussing ensemble 
(non-)equivalcnce in Sec. IV. 

III. MICROCANONICAL ENSEMBLE 

In the microcanonical ensemble, the natural starting 
point for thermodynamics is the entropy. For our pur- 
poses, we are interested in the thermodynamic limit of 
the configurational entropy, i.e., the contribution of the 
kinetic term in (2) to the entropy is disregarded. We will 
compute and discuss two different — but closely related — 
entropy functions, the first being a function of the poten- 
tial energy and the magnetization, 

s(v,m) = lim —In/ dtp 6 [v — Vw((p)] S [m — mjv(y)] , 

(6) 

where S denotes the Dirac distribution. The second en- 
tropy function we consider depends only on the potential 
energy, 

s(v) = lim -^-ln / dtpS [v - V N {(p)} , (7) 

N^oo iv J]R« 

and can be obtained from s(v, m) by maximization, 

s(v) = sup s(v, to). (8) 

m 

This leaves the calculation of s(v,m) as the main task 
of this section. Before doing so, we want to say a few 
words about how a phase transition manifests itself in 
the microcanonical entropy. 
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A. Microcanonical entropy and phase transitions 

In the canonical ensemble, a possible (and commonly 
accepted) definition of a phase transition is the presence 
of a nonanalyticity in the canonical free energy. In the 
microcanonical ensemble, conditions on the microcanon- 
ical entropy s for the existence of a phase transition are 
somewhat less established. An obvious way to start is 
to look for conditions on s which, in the case of ensem- 
ble equivalence, correspond to the definition of a phase 
transition from the canonical free energy density. Such 
conditions are 

1. the existence of an open subset of the domain of s 
on which the entropy is not strictly concave, or 

2. a nonanalyticity in s. 

Then it might appear reasonable to adopt the same con- 
ditions also for systems for which equivalence of ensem- 
bles does not hold. So in order to identify phase tran- 
sitions in the microcanonical ensemble, we will focus on 
concavity and analyticity properties of the entropy in the 
following. 



B. Large deviation theory and microcanonical 
entropy 

Large deviation theory is a branch of probability the- 
ory which is concerned with events of very low probabil- 
ity. This section provides an informal introduction as to 
how large deviation theory can be applied to calculate 
the microcanonical entropy. 

Let X N = (Xi, X2, . . • , Xn) be a sequence of N in- 
dependent and identically distributed random variables 
Xi G IR™ with empirical mean 



e N (x N ) 



X!+X 2 



x 



N 



N 



(9) 



and mean fi = E[JQ], where E[-] denotes the expectation 
value. The strong law of large numbers then states that 
the probability density to find &n =/= [i converges to zero 
as N goes to infinity. The theory of large deviations deals 
with the form of this convergence. Provided existence of 
the generating function E[e^* ,x ^], the rate function / is 
defined as the Lcgendrc-Fcnchcl transform 



I(x) 



sup 

/ 6R" 



(t,x) -lnE[e ( *' x) ; 



(10) 



where (■, •) denotes the Euclidean scalar product. Then 
Cramer's theorem (Theorem 1.4 in [20]) states that the 
probability P to find the empirical mean &n G [a;, i+da;] 
converges exponentially, 

P(e N G [x,x + dx]) = e- NI ^dx (11) 

in the limit N — ► 00, or, equivalently, 

-I(x)= lim — In P(e N G [x,x + dx]), (12) 



and 6^r is said to satisfy a large deviation principle 

The connection to statistical physics is made by ob- 
serving that the microcanonical entropy of a system with 
./V degrees of freedom, described by generalized coordi- 
nates Xi G R™, is defined as 



s N (x) 



1 

N 



In 



dx N s [ x -e N (x N )] . (13) 



The integral in this equation quantifies the volume 
in state space WL Nn occupied by microstates X N = 
(X\, X2, . ■ ■ , Xn) which are compatible with a certain 
macrostate x. If P(&n € [x,x + dx]) is the probability 
to find the system in a state with &n G [x,x + dx], the 
proportionality of this integral to P(&n & [x, x+dx]) fol- 
lows immediately from the assumption of equal a priori 
probability of the microstates, which is at the very basis 
of equilibrium statistical mechanics. The proportionality 
constant results in a physically irrelevant summand in 
the entropy and is omitted, leading to the expression 

s(x) = lim — In P(6 N G [x,x + dx]) (14) 

for the entropy in the limit of large N. Provided our 
problem satisfies a large deviation principle, the micro- 
canonical entropy s(x) for N — > 00 is identical to the 
negative rate function —I(x) from Eq. (10). Thus, the 
high-dimensional integral in (13) is reduced to one single 
integration [the expectation value E^*'^] in (10)] and 
a maximization, 



s(x) 



sup 

ten™ 



(t,x) -lnE[e< 4 <*>] 



(15) 



lnE[e<*' X> is strictly convex and infinitely many times 
differcntiable (sec proof of Theorem 1.4 in [20]). There- 
fore, there is at most one maximizer t x , and the derivative 
of the argument of the supremum in (15) has to be zero. 
As a consequence, t x is determined by the equation 



xM[e {t *' X) ] = E 



x) 



and the maximization can be rewritten as 



s(x) = -{t x ,x)+hiM[e^'^]. 



(16) 



(17) 



Note that —I(x), and therefore s(x), are smooth and 
strictly concave on their domain (Lemma 1.14 in [20]). 

In [21, 22], a more extensive discussion of the appli- 
cation of large deviation theory in statistical physics, 
including a detailed account of the properties of the 
rate function I{x), is given in a language well accessi- 
ble to readers with a physics background. A mathemat- 
ical treatment of large deviation theory can be found in 
[20, 23]. 



C. The if> -model without interaction 

We are interested in the entropy s(v,m) as a function 
of the potential energy v and the magnetization m as 
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defined in (6). As a first step, applying the concepts 
of large deviation theory, we compute a related entropy 
function s(z,m) in this section, where z is the variable 
associated with the on-site potential zn defined in (4). 
From s(z,m), we then obtain the desired s{v,m) by a 
simple transformation of variables. In the notation of 
Sec. IIIB, we choose x = (z,m) € R 2 . The adequate 
empirical mean Sat = (-zat^), tojv(</?)) is made up from 
(random) microscopic variables X{ = (\<p\ — <Pi) ■ 
We further assume uniformly distributed particles with 
probability density p{fi = ip) = Vy> G [— <p c , tp c ] and 
elsewhere. Then, the conditional equation (16) for t x = 
{t z ,t m ) leads to 



Z = \ Ui{t z ,t m ) - \ U2(t z ,t m ) 

to = wi(t z ,t m ) 
with parameter integrals 

/+£ c d^ fc c^-M^ 4 -^ 2 ) 



(18) 



(19) 



Inserting the solution t z {z,m), t m (z,m) of the system 
of equations (18) into (17) yields the microcanonical en- 
tropy for the noninteracting <p 4 -model, 

s(z, to) = — t z (z, m)z — t m (z, m)m 

+ ln / dipe tmV+t *(^ 4 -^ 2 ), (20) 



where we let ip c tend to infinity in the following. 

An analytic solution of (18), and therefore of the en- 
tropy s in (20), is not possible in general. Instead, we will 
state — without proof — some of the properties of s. From 
the definition (4) of the on-site potential one obtains im- 
mediately z G [—4, oo) for the range of z. For any fixed 
z > -\, values in the range [-m max (z), +TO max (z)] are 
accessible for the magnetization to. The extremal values 
nimax(z) are defined by the equation 



1 4 1 2 

Z = — 711 — — 77? 

A max o max 



(21) 



and correspond to states in which all degrees of freedom 
ifi take on the same value. This characterizes the domain 
of s as illustrated in Fig. 1. 

According to the previous section, s(z,m) is strictly 
concave and analytic on its domain. Taking into account 
the symmetry of the model [invariance of zn(p) under 
substitution ip — > —(f], it follows that s(z,m) is maxi- 
mal at m = for any given z ^ — j (see Fig. 1 for an 
illustration) . 



D. Microcanonical entropy of the mean- field 
y 4 -model 



So far we have considered the mean-field ip -model 
without intcrparticle interaction. Turning on this inter- 
action changes the potential energy of the system from 




FIG. 1: Contour plot of the entropy s(z,m) on its domain 
(gray) from a numerical evaluation of (20). The entropy is a 
concave function and minimal on its boundary for maximum 



Note that the bold line at z — 



is not a contour. 



Zn to Vn = Zn 



2 " L N- 



This relation allows to derive 



the desired entropy function 

s(v, to) = s (v + 



TO 2 , TO,) 



(22) 



of the interacting system from the above computed 
s(z, to) by a simple transformation of variables. 

The smooth transformation of variables in (22) retains 
the analyticity of s(z, to), so s(v, to) is again an analytic 
function. In order to sec whether the concavity prop- 
erty is conserved as well, we distinguish the two cases of, 
respectively, negative and positive coupling J. 

J < 0: Let ,fi/ denote the Hessian of some function 
/. Using (22), we can calculate the determinant 



detf) s = detSjg + J 



d 2 s(z, to) ds(z, to) 



dz 2 



d z 



(23) 



< 



The concavity of s implies that det > and — 
0. If we further take into account that ds ^- m ^ = —f > o 

o z 

it follows that dct.f) s > for J < 0, therefore Sj s is 
negative definite and s(v,m) is again a strictly concave 
function. From the symmetry of the model it follows 
that, for any fixed value of v, there is a unique maximum 
at to = 0. Then the entropy as a function of v only, 
s(v) = sup m s(w,TO) = s(v, 0), is likewise analytic and 
concave, and no phase transition is found to take place 
for anti-ferromagnetic coupling J < 0. In the following 
we will focus on the more interesting ferromagnetic case. 

J > 0: For arbitrary positive coupling J, the domain 
of s is no more a convex set (see Fig. 2), so s is not a 
concave function, and we expect a phase transition to 
take place. 

Numerically, we find that for fixed values of v above a 
critical value v Cl s(w,to) attains its maximum at magne- 
tization to = 0, whereas below v c there is a minimum at 
to = and two maxima occur at nonvanishing values of 
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FIG. 2: Contour plot of the entropy s(v,m) for coupling con- 
stant J = 1. The gray (light and dark) hatched area is the 
domain of s, dark gray is used for what is obtained from a 
deformation (variable transform) of the part visible in Fig. 1. 
The dashed line marks the position of the maximum of s with 
respect to m for every fixed value of v. 



to (Fig. 2). In order to track down the transition point, 
we expand the entropy s in to around to = 0, yielding 



s(v, to) = s(v, 0) 
- \ (jr z (v 1 0) + 



1 



w 2 (r>,0),0) 



to 2 + C(to 4 ), (24) 



where T z (v,m) = t z {v + ^m 2 ,m). Here, O denotes Lan- 
dau's order symbol, and the integral 



i-i(-^)+i f (-^r 



(25) 



as defined in (19) can be expressed in terms of modified 
Bessel functions of the first kind Ik ■ The transition point 
is marked by a change of sign of the coefficient of the 
quadratic term in (24), so the critical energy is given 
via Eq. (18) by v c = z(t Zc ,0), where t Zc is determined 
implicitly by 



Jt Zc ui 2 (Tz o ,0) = -1. 
Making use of this equation and the identity 



wa(*.,0)-i 



the critical energy can be expressed as 
v c (J) 



1_1 

4r, c (J) 



(26) 



(27) 



(28) 



where the dependence of v c and t Zc on the coupling con- 
stant J has been indicated. We can immediately read 
off that v c (l) = and, since t Zc < 0, that v c > for 
all J > 1. Solving explicitly for t Zc or v c appears to be 
infcasible, but we can derive asymptotic expansions for 
large and small coupling J, respectively 



Small J > 0: The asymptotic expansion of (25) for 
large — r 2e yields 



^2{T Zc ,0) 



-4t, 



-4r z 



(29) 



Inserting this result into (26) and (28), an expansion of 
the critical energy in the limit of small coupling J can be 
obtained. 



vc(j) = -! + ! j + o(j 2 ). 



(30) 



Large J > 0: Analogously to the preceding para- 
graph, expansion of (26) for small negative t Zc yields an 
expression for the critical energy in the limit of large J, 



{J)=a 2 J 2 -{2a 2 -\)J- 



4 



1 + 

8 ' 64 a 2 



(31) 

with a = r(|)/r(i) w 0.338. 

In Fig. 3, the numerically computed behavior of the 
critical energy v c as a function of the coupling J is con- 
fronted with the asymptotic expansions (30) and (31). 
Note that, for large enough J, the critical energy takes 




-0.4 i 1 1 1 1 1 1 1 H J 

0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

FIG. 3: Critical energy v c as a function of the coupling con- 
stant J from numerical computation (dots), confronted with 
the two asymptotic expansions (30) and (31). 

on positive values, as in fact it is divergent in the limit 
of large J. This observation will be of importance for 
the discussion of the relation between phase transitions 
and the topology of configuration space submanifolds in 
Sec. V. 



E. Reduced entropy s(v) 

In Sec. HID, the entropy s(v,m) was found to be 
an analytic function. If we release the control param- 
eter m, the system will maximize its entropy to s(v) = 
sup m s(v,m), and we will discuss the analyticity proper- 
ties of this latter function in the present section. Since 
s(v, to) is minimal at its boundary, a vanishing first 
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derivative with respect to m is a necessary and suffi- 
cient condition for the extremalization. From (20) and 
(22) we find the extremizing values of the magnetization 
rh = LUi(T z ,r m ) to be determined by 

wi(f z ,f m ) = — (32) 

J T z 

The numerical solution of this equation is plotted in 
Fig. 4. As expected from symmetry considerations, 

7~m 




Tz 



-1.6 -1.4 -1.2 -1.0 -0.8 



FIG. 4: Numerical solution of the extremalization condition 
(32) for J = 1. f m = is a solution for any f z < 0. Below 
the critical value r Zc there are two further solutions. 

Eq. (32) is trivially solved by f m = Vf 2 £ 1R _ , so rh = 
is always an extrcmum. The Taylor expansion in (24) re- 
vealed that this solution corresponds to a maximum for 
t z > t Zc and to a minimum for f z < r Zc . Numerics as- 
certain that, for f z > t Zc , the maximum at rh = is 
a global one, and the system is found to be in a para- 
magnetic phase. For f z < r Zc two additional solutions 
±T m (T z ) appear. These solutions correspond to two max- 
ima at rh = — ^J'fr ^ 7^ and reveal the system to be in 
a ferromagnetic phase. The transition between the two 
phases is a continuous one (see dashed line in Fig. 2). The 
critical exponent (3 governing the asymptotic behavior of 
the magnetization in the vicinity of the transition point 
can be calculated analytically and is found to be j3 — |. 
In the Appendix, we show that the occurrence of clas- 
sical (=mean-field) values of the critical exponents can 
be traced back to the analyticity of the entropy function 
s(v,m), not only in the case of the mean-field <y9 4 -model, 
but generically for mean field-like systems. 

In conclusion, the entropy s(v) as a function of the 
potential energy only is given by 

v n fs(v,0) for v > v c , 

S ( V > = 1 ( r m (T,H , ^ (33) 

[S(V, J^fT 1 ) t0TV<V c , 

where r z is determined by 

v = iw 4 (T z ,r m (r z )) - iw 2 (r z ,T m (r z )) - . (34) 

A nonanalyticity corresponding to a continuous phase 
transition occurs at v = v c . 



IV. ENSEMBLE (NON-)EQUIVALENCE 

Results for the mean-field </j 4 -modcl from computa- 
tions within the canonical ensemble have been reported 
in the literature, showing the system to undergo a contin- 
uous phase transition [16, 18, 19]. Under suitable condi- 
tions, the different Gibbs ensembles, like the microcanon- 
ical and the canonical one, are known to lead to identical 
numerical values for the typical system observables of 
interest in the thermodynamic limit N — > oo. In this 
case, one speaks of equivalence of ensembles. For a large 
class of systems with short-ranged interactions, equiva- 
lence of ensembles is known to hold [24, 25]. A system 
with mean field-like interactions like the one considered 
here is clearly not in this class. Nevertheless, compar- 
ing the microcanonical result for the critical energy v c 
obtained in the previous section with the canonical one 
reported in [16, 19], we observe that the expressions are 
identical. This is not at all unexpected, as the theorem 
of Gartner and Ellis [23] guarantees equivalence of en- 
sembles (on the thermodynamic level we are interested 
in here, see [26] for a precise definition) whenever there 
is no discontinuous phase transition in the canonical en- 
semble [22]. 

Despite this ensemble equivalence, the information en- 
coded in the microcanonical results is richer than that 
of canonical results: The convex region of s (the interior 
of the dashed parabola-like curve in Fig. 2) is inacces- 
sible within the canonical ensemble, even in the pres- 
ence of a magnetic field. As a consequence, the analyt- 
icity of the microcanonical entropy s(v,m) — which will 
be essential for our discussion of the relation between 
phase transitions and the topology of configuration space 
submanifolds can not be inferred from canonical results. 



V. TOPOLOGY OF CONFIGURATION SPACE 
SUBMANIFOLDS 

The topology of the configuration space submanifolds 
M v as defined in (1) has been analyzed for the mean- 
field </> 4 -model in Refs. [14-16] in the framework of Morse 
theory. To this aim, critical points [32] of the potential 
energy function (3) were determined, i.e., points tp c £ 
1R, N at which the exact differential of vn vanishes, 

dv N (pc) = 0. (35) 

For a given coupling constant J, all critical values vn{^c) 
were found to lie within a finite interval, 

M<Pc) € [-|(J+1) 2 ,0] , (36) 

independently of the number N of degrees of freedom. 
The observation uat(</J c ) ^ is already sufficient for our 
purposes, as the non-critical neck theorem [27] then guar- 
antees that the topology of the M v remains unchanged 
for all v > 0. Recalling the observation from Sec. HID 
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that the critical potential energy v c of the phase tran- 
sition can attain arbitrarily large values, well above the 
zero upper bound of the critical values of vn, we are led 
to conclude that the phase transition is found to be un- 
corrected with any of the topology changes occurring in 
the family {M v } of configuration space submanifolds. 

VI. SINGULARITY-GENERATING 
MECHANISMS 

The reason why, as observed in the preceding section, 
the phase transition of the mean-field <yS-model need not 
be related to the topology changes in the M v is the 
existence of a further singularity-generating mechanism. 
This mechanism is the maximization over the magneti- 
zation m in Eq. (8) of the real analytic entropy func- 
tion s(v,m), resulting in a nonanalytic s(v). In the 
present section, we will explain why this mechanism is 
restricted to systems with long-range interactions. Fur- 
thermore, it is argued that a topology change in the M v 
generates a thermodynamic singularity on a more funda- 
mental level. Although a heuristic argument, we believe 
that it illustrates the essential difference between the two 
singularity-generating mechanisms discussed. For the 
moment, we will restrict ourselves to continuous phase 
transitions, but modifying the argument as to apply to 
discontinuous transitions is straightforward. 

In Fig. 2, the typical shape of the entropy s(v, m) of a 
long-range system in the presence of a symmetry break- 
ing continuous phase transition with some order parame- 
ter m is shown: two maxima in s(v, m) with respect to m 
at potential energies v below the critical energy are op- 
posed to a single maximum at potential energies above 
the critical energy. As illustrated in Sec. Ill for the mean- 
field </3 4 -model, even an analytic entropy function s(v, m) 
can entail such a phase transition. 

This is different for a system with short-range interac- 
tions [33]. In this case, the entropy as a function of some 
macroscopic variables v, m, x,... is known to be a con- 
cave function [24], and the typical shape of s(v, m) in the 
presence of a phase transition is obtained by taking the 
concave hull of a long-range system's entropy (see Fig. 5 
for an illustration). Due to its concavity, the entropy s 
is nonanalytic not only as a function of v only, but the 
same holds true for s(v, m) or for s(v, x) for some typical 
macroscopic quantity x. 

We will try to trace back the difference between the 
long-range and the short-range case to the integrals gov- 
erning the entropy functions. We use a co-area formula 
from [28] to rewrite Eq. (7) as 



with norm 
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FIG. 5: Sketch of a typical entropy function s(v, m) of a sys- 
tem with short-range interactions showing a continuous phase 
transition. The shape is qualitatively that of the concave hull 
of the entropy of a long-range system (like the one depicted 
in Fig. 2). The dashed line marks the spontaneous magneti- 
zation. 

dM v = E„ denotes the boundary of M v . Similarly, the 
entropy as a function of v and x [with x = m yielding 
Eq. (6) as a special case] can be expressed as 

s(v,x)=lim— — — j| S [x - x N (ip)) , (39) 

N^oo N J dMv \\VV N \\ 

with some smooth function x n ■ K." — ► M corresponding 
to the macroscopic variable x. Then the above discussed 
analyticity properties of the entropy functions have the 
following implications for the integrals in (37) and (39): 

For the case of short-range interactions, we observed 
above that, in the presence of a phase transition, s(v ) as 
well as s(v,x) are necessarily nonanalytic, regardless of 
the particular form of xpj. So if its not the form of the 
integrand which is crucial, we are led to conjecture that 
it is the domain of integration dM v , and in particular 
some topology change in the M v , which accounts for the 
nonanalyticities in s(v) and s(v,x). This argument is 
consistent with the theorem in Refs. [12, 13], proving the 
necessity of a topology change for a phase transition to 
take place for a certain class of short-range systems. 

For systems with long-range interactions, we found 
that a nonanalytic s(v) can emanate from an analytic 
s(v,m). In this case, as opposed to our reasoning for 
short-ranged systems, the emergence of a nonanalyticity 
in the thermodynamic limit TV — > oo appears to depend 
on the shape of the integrand, not on the domain of inte- 
gration dM v . This explains on a heuristic level why, for 
systems with long-range interactions, a phase transition 
is not necessarily connected to a topology change in the 
family {M v } of configuration space submanifolds. 

The major difference between the two singularity- 
generating mechanisms is that a topology change in {M v } 
may lead to a nonanalyticity, regardless of the particular 
choice of arguments of the entropy function, which is ob- 
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viously not true for a singularity stemming from a maxi- 
mization over one of the variables of an analytic entropy 
function. It is in this sense that by a topology change a 
nonanalyticity is generated on a more fundamental level. 



VII. SUMMARY AND CONCLUSIONS 

Applying a large deviation technique, we have obtained 
an exact expression for the microcanonical (configura- 
tional) entropy s of the mean-field (p 4 -model as a func- 
tion of the potential energy v and the magnetization m. 
Although the system undergoes a continuous phase tran- 
sition at some critical energy v c , the entropy s(v,m) is 
found to be an analytic function in both variables. Only 
the reduction s(v) = sup m s(v, to) to an entropy function 
of one variable gives rise to a nonanalyticity, correspond- 
ing to the phase transition of the model. 

As expected, the phase transition is found to be gov- 
erned by classical (=mean-field) critical exponents a = 
0, (i = |, . . . The occurrence of these values can be traced 
back to the analyticity of the entropy s(v, to), not only in 
the case of the mean-field </? 4 -model, but also generically 
for systems with mean field-like interactions. 

A topology change within the family {M v } of configu- 
ration space submanifolds M v [as defined in (1)] is known 
to be a mechanism which can give rise to a nonanalyt- 
icity in some thermodynamic potential, thus entailing a 
phase transition. The observed analyticity properties of 
the entropy functions allow us to individuate the max- 
imization over one variable of an analytic function as a 
further such mechanism. For the mean-field <y9 4 -modcl, 
the topology changes were found to be unrelated to the 
phase transition. This, previously unexpected, obser- 
vation can be explained by the presence of this second 
singularity-generating mechanism. Concavity of the en- 
tropy of short-range systems restricts the mechanism of 
singularity generation by maximization to systems with 
long-range interactions. 

We want to conclude with a remark on the implica- 
tions the results of this article have for what we called 
the topological approach to phase transitions in the In- 
troduction. The identification of a topology change in 
the M v as one out of two, or several, mechanisms to 
produce a thermodynamic singularity alters profoundly 
the perspective from which the topological approach has 
to be viewed. Clearly, in its generality, the above men- 
tioned topological hypothesis [2] has to be discarded: in 
general, a topology change in the M v at some potential 
energy v = v c is not necessary for a phase transition to 
take place at v c . But this docs not diminish the inter- 
est in the topological approach. Instead, the classifica- 
tion of phase transitions according to the mechanism by 
which the thermodynamic singularity is generated opens 
up as a perspective, with topology changes in configu- 
ration space being one — for short-range systems possi- 
bly the most prominent one — among these mechanisms. 
In addition to the two singularity-generating mechanism 



discussed in this article, there are indications for at least 
one further such mechanism: also in short-range systems 
with non-confining potentials, a topology change does 
not appear to be necessary for the existence of a phase 
transition [29], while at the same time the maximization 
mechanism is excluded due to the short-rangedness. One 
might conjecture that it is the fact that, as a consequence 
of the potentials being non-confining, the domains of in- 
tegration in Eqs. (37) and (39) are non-compact mani- 
folds, which might render a third singularity-generating 
mechanism possible, but this point requires further in- 
vestigation. 
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APPENDIX: CRITICAL EXPONENTS FOR 
MEAN FIELD-LIKE SYSTEMS 

In this appendix, we show that generically (in a 
sense which will become clear in the following) classical 
(=mean-field) values follow from the analyticity of the 
entropy s. The mean-field (p -model can then be proved 
explicitly to be such a generic case by showing that cer- 
tain coefficients of a Taylor expansion of s do not vanish. 

We begin with the remark that the arguments of 
Sec. Ill, guaranteeing the analyticity of the entropy 
s(v,m), apply not only to the mean-field y? 4 -model, but 
to general mean field-like systems for which it is possible 
to express the potential energy by means of macroscopic 
variables of the kind X)i=i 9k{fi) with some functions 
[34]. In fact, analyticity of the entropy holds even for 
systems with more realistic long-range forces, whenever 
the potential energy can be expressed by macroscopic 
variables plus some remainder, vanishing in the thermo- 
dynamic limit (see [21] for an example). 

In order to simplify the presentation, the following 
derivation is formulated for systems which are symmet- 
ric with respect to a (scalar) order parameter to, i.e., 
s(v, to) = s(v, —to), but a more general setting is possi- 
ble. Assuming analyticity, s(v,m) can be expanded into 
a Taylor series around to = 0, 

s(v,m) = s{v,Q) + f 2 (v)m 2 + f 4 {v)m i + 0{m e ), (A.l) 
with Mv) = #4r| n and Uv) = #4r| n . Odd or- 

I dm 2 I m=0 1 am 4 1 777=0 

ders in m vanish due to the system's symmetry. We will 
understand by the above mentioned generic case that no 
coefficients in the expansion (with respect to both, to and 
v) vanish accidentally, so we assume nonzero coefficients 
throughout. 

The cxtrcmizing values of the magnetization are deter- 
mined by gjjj = 0. Applying this to (A.l) leads to the 
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trivial solution rh — and, for small fh =/= 0, to 

™(v) ~ iv^M (A.2) 

whenever the radical is positive. For the concept of criti- 
cal exponents to be meaningful, we assume a continuous 
phase transition to take place. This implies /2(f) to be 
zero at some critical point v = v c , and the expansions of 
$1 and fi around v c are of the form 

h{v) = &(v c ) (v - v v ) + O ((v - v c ) 2 ) (A.3) 
h{v) = U(v c ) + 0(v-v c ) (A A) 

Inserting these expressions into (A.2), we obtain 

m(v) ~ ±^J-^^y/v -v c <x(v- v c ) 1/2 (A.5) 

in leading order for the equilibrium magnetization in the 
vicinity of the v c . Similarly, the asymptotic behavior of 



other thermodynamic quantities can be determined from 
a generic analytic entropy function (see [30] for deriva- 
tions). To obtain critical exponents, the proportionali- 
ties to v — v c have to be translated into proportionalities 
to the reduced temperature t = T ~ijF c by means of the 
asymptotic relation v — v c ~ t 1 ~ a , where a is the criti- 
cal exponent of the specific heat. In this way, we finally 
obtain a critical exponent (3 = i for the equilibrium mag- 
netization of a system with analytic entropy s undergoing 
a continuous phase transition. 



An elaborate but straightforward calculation confirms 
that, for the mean field ip 4 -mode\, the above conditions 
(A.3) and (A. 4) on f 2 and / 4 hold with f^fi ^ 0, ren- 
dering this model a representative of the class of generic 
systems with classical critical exponents. 
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